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Abstract. We consider effect of stochastic sources upon self-organization process 
being initiated with creation of the limit cycle. General expressions obtained are 
applied to the stochastic Lorenz system to show that offset from equilibrium steady 
state can destroy the limit cycle at certain relation between characteristic scales of 
temporal variation of principle variables. Noise induced resonance related to the limit 
cycle is found analytically to appear in non-equilibrium steady state system if the 
fastest variations displays a principle variable, which is coupled with two different 
degrees of freedom or more. 
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1. Introduction 

Interplay between noise and non-linearity of dynamical systems [1] is known to arrive at 
crucial changing in behavior of systems displaying noise-induced [21 [3] and recurrence 
[H [5] phase transitions, stochastic resonance P, [7], noise induced pattern formation 
[HI [9], noise induced transport [IHl E] et cetera (see Ref. [H], for review). The 
constructive role of noise on dynamical systems includes hopping between multiple 
stable attractors [T2l[T3] and stabilization of the Lorenz attractor near the threshold of its 
formation [HI [15] . Such type behavior is inherent in finite systems which involve discrete 
entities (for instance, in ecological systems individuals form population stochastically in 
accordance with random births and deaths). Examples of substantial alteration of finite 
systems under effect of intrinsic noises give epidemics [I6]-[T8], predator-prey population 
dynamics [19], [20], opinion dynamics [21], biochemical clocks [22|,[23], genetic networks 
[21], cyclic trapping reactions [25], etc. 

Within the phase-plane language, phase transitions pointed out present the simplest 
case, where a fixed point appears only. We are interested in studying more complicated 
situation, when the system under consideration may display oscillatory behavior related 
to the limit cycle appearing as a result of the Hopf bifurcation [261 EZ]- It has long 
been conjectured [28] that in some situations the influence of noise would be sufficient 
to produce cyclic behavior [29]. Recent consideration [30] allows the relation between 
the stochastic oscillations in the fixed point phase and the oscillations in the limit cycle 
phase to be elucidated. Moreover, excitable [3l], bistable [32] and close to bifurcations 
[33] systems display oscillation behavior, whose adjacency to ideally periodic signal 
depends resonantly on the noise intensity [31] , that was a reason to call these oscillations 
coherence resonance |3l] or stochastic coherence [TT]. Control of the coherence resonance 
regime was shown to be achieved with a time-delayed feedback, which enables to increase 
or decrease the regularity of motion [35]. Characteristically, a quasioscillatory behavior 
may be organized without any input signal, provided a stochastic nonlinear system has 
itself an intrinsic time scale. If this scale is driven by a multiplicative noise, which induce 
bistable behavior in a deterministically monostable medium, then a doubly stochastic 
resonance arises [36] . 

The simplest way to formulate the model related to systems with finite number 
< oo of constituents is to consider the sum S = Y^iLi 6t oi random state vectors 
with components a = 1, . . . ,d. Then, the state vector 

S = NX + VNx (1) 

is decomposed into a deterministic component being proportional to total system size 
and a random one to be proportional to its square root [37]. In the limit of infinite 
particle numbers ^ oo, such systems are faithfully described by deterministic 
equations to find time dependence X(t), which addresses the behavior of the system 
on a mean-field level. On the other hand, a systematic study of corrections due to 
finite system size can capture the behavior of fiuctuations x{t) about the mean-field 
solution. These fiuctuations are governed with the Langevin equations, however, in 
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difference of approach |30j, we consider multiphcative noises instead of additive ones, 
on the one hand, and nonhnear forces instead of hnear ones, on the other. Within 
such framework, the aim of the present paper is to extend analytical descriptions [30] of 
finite-size stochastic effects to non-equilibrium systems where noises play a crucial role 
with respect to periodic limit cycle solution creation or its supression. We will show that 
character of the stationary behavior of non-equilibrium system is determined by relation 
between scales of temporal variation of principle variables as well as their coupling. In 
contrast to the doubly stochastic resonance [36], we consider the case when multistable 
state is caused by both multiplicative noise and offset from equilibrium state. 

The paper is organized along the following lines. In Section [2], we obtain conditions 
of the limit cycle creation using pair of stochastic equations with nonlinear forces and 
multiplicative noises. Sections [HI H] are devoted to consideration of these conditions 
on the basis of stochastic Lorenz system with different regimes of principle variables 
slaving. According to Section [3] the limit cycle is created only in the case if the most 
fast variation displays a principle variable, which is coupled nonlinearly with two other 
degrees of freedom or more. Opposite case is studied in Section H] to show that the limit 
cycle disappears in non-equilibrium steady state. Section [5] concludes our consideration. 



2. Statistical picture of limit cycle 

According to the theorem of central manifold |26j, to achieve a closed description of a 
limit cycle it is enough to use only two variables Xa, a = 1,2. In such a case, stochastic 
evolution of the system under investigation is defined by the Langevin equations [39] 

Xa = &^+GaCa{t), a =1,2 (2) 

with forces f^°''^ = f^°'\xi,X2) and noise amplitudes Qa = Ga{xi,X2), being functions of 
stochastic variables Xa, a = 1,2; white noises Ca{t) are determined by usual conditions 
{(ait)) = 0, (Ca(^)C/3(^')) = ^a/s^it — t')- Within the assumption that microscopic transfer 
rates are non correlated for different variables 

ex. ^ SGG below) , the probability distribution 
function V = V{xi,X2] t) is determined by the Fokker-Planck equation 

where components of the probability current take the form 

jia) ^ _ 1 ^ ^ {gaQf^V) (4) 

2 ^^-^ C)Xf3 

with the generalized forces 

^H = /(°) + AV^%^, (5) 
^1 dxp 

being determined with choice of the calculus parameter A G [0, 1] (for Ito and 
Stratonovich cases, one has A = and A = 1/2, respectively). Within the steady state. 
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the components of the probabihty current take constant values jjf^ and the system 
behaviour is defined by the following equations: 

^ (^'^) + ^ ^^^^^^^ - 2-^^'^^ = -2^°'^' 

^ iQig2V) + ^ {g'^p) - 2p'^v = -2jf . (7) 

Multiplying the first of these equations by factor ^2 and the second one by ^1 and 
then subtracting results, we arrive at the explicit form of the probability distribution 
function as follows: 

D(Xi,X2) = (^2^(^)- ^1-^(2)) (8) 



1 

+ 2 



^1 n ~ ^1^2 



dxi 8x2 ) \ dxi 8x2 

This function diverges at condition 



that physically means appearance of a domain of forbidden values of stochastic variables 
Xai which is bonded with a closed line of the limit cycle. Characteristically, such a line 
appears only if the denominator T'(xi, X2) of fraction ([8]) includes even powers of both 
variables xi and X2lil 

It is worth to note that the analytical expression ([8]) of the probability distribution 
function becomes possible due to the special form of the probability current (jlj), where 
effective diffusion coefficient takes the multiplicative form Dap = gagp- In general case, 
this coefficient is known to be defined with the expression [ID] 

Vap = Y.iai>glg% (10) 

ab 

where kernel lab determines transfer rate between microscopic states a and 6, whereas 
factors Qa and (7^^ are specific noise amphtudes of values Xa related to these states. We 
have considered above the simplest case, when the transfer rate lab = / is constant for 
all microscopic states. As a result, the diffusion coefficient (fTO!) takes the needed form 
T^af3 = gag 13 with cumulativc noise amplitudes = \/lY,a9a ^^"^ = ^HbQ^p- 

3. Noise induced resonance within Lorenz system 

As the simplest and most popular example of the self-organization induced by the 
Hopf bifurcation, we consider modulation regime of spontaneous laser radiation, whose 
behaviour is presented in terms of the radiation strength E, the matter polarization P 

\ Archetype of closed curves presents the circle x\ + x\ — 1. 
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and the difference of level populations 5" ^7\. With accounting for stochastic sources 
related, the self-organization process of this system is described by the Lorenz equations 

teE = [-E + geP - ^{E)] + QECit), 

TpP={-P + apES)+gpat), (11) 

TsS=[{S,-S)-asEP]+gsm. 

Here, overdot denotes differentiation over time t; te^p^s and aE,p,s > are time scales 
and feedback constants of related variables, respectively; gE,p,s are corresponding noise 
amplitudes, and Se is driven force. In the absence of noises {qe = Qp = Qs = ^) and at 
relations tp,ts ^ te between time scales, the system (fTT!) addresses to limit cycle only 
in the presence of the nonlinear force [H] 



characterized with parameters k > and E^. In this Section, we consider noise effect 
in the case of opposite relations te <^ Tp, ts of time scales, when periodic variation of 
stochastic variables becomes possible even at suppression of the force (fT2!) . 

It is convenient further to pass to dimensionless variables t, C, E^ P, S, qe, gp, gs 
with making use of the related scales: 

Tp; Cs = Tp^^"^; Es = (apas)"^/^ Ps = (alapas)'^^"^, Ss = (a£;ap)~^; 

9e = {rp/apasf^, gp = {tp/ a^apasf^ , g's = Tp^'^/aEap. (13) 

Then, equations ( ITTl) take the simple forn]§| 

a-^E = -E + P-ip{E)+gEat), 

P = -P + ES + gpat), (14) 

ie/a)S=iSe-S)-EP + gsat), 

where the time scale ratios 

(t = tp/te, e = Ts/TE (15) 

are introduced. In the absence of the noises, the Lorenz system f|T^ is known to show 
the usual bifurcation in the point Se = 1 and the Hopf bifurcation at the driven force 
[381 [371 

= ^ (16) 

However, the noiseless limit cycle {gE = gp = gs = ^) is unstable and the Hopf 
bifurcation arrives at the strange attractor only. 

With switching on the noises, the condition te <^ Tp allows for to put l.h.s. of the 
first equation ( iMj) to be equal zero. Then, the radiation strength is expressed with the 
equality 

E = P + gEat),m (17) 

§ These equations are reduced to the initial Lorenz form 38^ if we set X = \J a jeE, Y = ^cr/eP, 
Z = Se - S, r = Se, b= a/e and = gp = gs ^ 0. 



Noise induced oscillations in non- equilibrium steady state systems 



6 




Figure 1. Steady state distribution function 
Se = 0.5, gE = 0.5, gp = 1.376, gs = 2.5 



at J, 



iP) 



1, ^0 



10, Tp = TS, 



whose insertion into the system flMj) reduces it into two-dimensional form 

p = -Pii-s) + gpat), 

s = ia/E) [iSe-s)-p']+gsat) 

with the effective amphtudes of muhiphcative noises 

Gp = ^Jgl + 9lS\ Qs={rp/Ts)^jgl + glP^ 
and the generahzed forces 



{9s/9E? + P^ 



rs/rp \{gp/gEf + S^- 
jr(s) = (^rp/rs) [(Se - S) 
In this way, the probabihty density ([H]) takes infinite values at condition 



P' 



{9p/9e 


)2 + S2 


{9s/9e) 


2 + P2 



9l<y ( 91 



9l 



9l 



^ + S^P{1 -S)+\^- + S 



\9l 



2 £ V^l y 



2 \9e / 



9l 



(18) 



(19) 



(20) 



ii|L + P2 (S-, - S) - p2 
\ 9e 



P = 0, 



(21) 



where we choose the simplest case of the Ito calculus (A = 0). 

Reduced Lorenz system f|T8|) has two-dimensional form ([2]), where the role of 
variables xi and X2 play the matter polarization P and the difference of level populations 
S. According to the distribution function ([8]) shown in FigUl the stochastic variables 
P and 5* are realized with non-zero probabilities out off the limit cycle only, whereas in 
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Figure 2. Form of limit cycle determined with Ea. (PT|) a.t e — 1, a — 1 and: 

a) gE — 0.5, gp — 11, gs — Q (curves 1-3 relate to Se = 0.5,1.0,2.0, respectively); 

b) Se = 0.5, gp = 7.5, gs — 6.5 (curves 1-3 relate to gs = 1.0,0.6,0.5, respectively) 




5 10 15 5 10 15 Qp 



Figure 3. Phase diagrams of the limit cycle creation at £ = 1, cr = 1 and: a) gs = 0.5, 
curves 1-4 correspond to Se = 0.0,0.5,1.0,2.0, respectively; b) Se = 0.5, curves 1-3 
correspond to gs = 0.1, 0.5, 1.0, respectively (diamonds relate to the values gp and gs, 
for which limit cycles in Figl2] are depicted) 

its interior the domain of forbidden values P, S appears. That is principle difference 
from the deterministic limit cycle, which bounds a domain of unstable values of related 
variables. The form of this domain is shown in FigI2] at different values of the noise 
amplitudes qe, gp, gs and driven force S^. It is seen, this domain grows with increase of 
the driven force Se, whereas increase of the force fluctuations gE shrinks it. On the other 
hand, phase diagrams depicted in Figj3] show that increasing the noise amplitudes of 
both polarization and difference of level populations enlarges domain of the limit cycle 
creation (more exactly, the noise amplitude gE shrinks this domain from both above 
and below, whereas increase of the driven force makes the same from above only). 

The principle peculiarity of the limit cycles obtained is that their form, determined 
with Eq. (l2T]) . does not depend on a non-equilibrium degree fixed by the stationary 
probability currents Jq , whereas the probability ([8]) itself does not equal zero at 
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conditions Jq ^ only. In this connection, one should be pointed out a non- 
triviality of the problem of numerical solution of the reduced Lorenz system ( ITSl) . 
which determines these limit cycles initially. Indeed, resolving this problem proposes 
the following steps: i) direct solution of the stochastic equations (fTSll to find a set of the 
time dependencies P{t) and S{t)] ii) numerical determination of the time-dependent 
probability V{P,S;t) to realize entire set of possible solutions of the equations ( TT8l) : 
iii) selection of non-equilibrium solutions, which obey to the steady state condition 
= Jq°'\ a = 1,2, determined with the probabillity current @; iv) calculation of 
the probability distribution Vs{P, S) of the steady state solutions; v) determination of 
the stochastic limit cycle according to the condition Vs{P, S) = oo. Realization of this 
program is in progress. 



4. Lorenz system without limit cycle 



According to Ref. [H], at conditions Tp -C te, ts, the deterministic system {gE,p,s = 0) 
has a limit cycle only at large intensity k of non linear force (fT2l) . In this case, it is 
convenient to measure the time t in the scale te and replace Tp by te in set of scales 
( fT3l) . Then, one obtains instead of Eq.( fT7|) the relation 

P = ES + gpat), (22) 
due to which the Lorenz system f|T^ is reduced to two-dimensional form 

E = -[E{l-S) + v{E)]+gEat), 



s 



s,-s{i + E') +gsat) 



with the effective noise amplitudes 

Ge = \/9p + 9e, = e 
The generalized forces are as follows: 

T^E) = -[E{l-S) + ^{E)], 



9l 



glE\ 



(23) 



(24) 



+ X^E. 



1 + (gE/gp) 



e \{gs/gpf + E^- 



(25) 



The probability distribution function ([8]) diverges at condition 

\2 , T^2 



{gs/gpY + E ^ 
1 + (gE/gp) 



ME) + E{1-S)] 



{gs/gpY + E^ 



Se - 3(1 + E^ 



2 



g'pE 



0, 



(26) 



^ l + iQE/gp) 

being the equation, which does not include even powers of the variable S. 

As a result, one can conclude that offset from equilibrium steady state destroys 
a deterministic limit cycle at the relations Tp -C te, ts between characteristic scales. 
This conclusion is confirmed with FigJH which shows divergence of the probability 
distribution function on the limit cycle of variation of the radiation strength E and 
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Figure 4. Steady state probability distribution function in dependence of the radiation 
strength E and the difference of level populations S at conditions Tp te ^ ts, 
K = 10, Se = 11.6, gE = 0.2, gp = 0.2, gs = 0.2 and different probability currents 
Jq'^-* , Jq'^-* (shown on the panels related) 

the difference of level populations S at zeros probability currents Jq and Jq only. 
With increase of these currents the system escapes from equihbrium steady state and 
maximum of the distribution function shifts to non-closed curves to be determined with 
equation (J26il . 

5. Conclusion 

We have considered effect of stochastic sources upon self-organization process being 
initiated with creation of the limit cycle. In Sections [3l HI we have applied general 
relations obtained in Section [2] to the stochastic Lorenz system. We have shown that 
offset from equilibrium steady state can destroy or create the limit cycle in dependence 
of relation between characteristic scales of temporal variation of principle variables. 

Investigation of the Lorenz system with different regimes of principle variables 
slaving shows that additive noises can take multiplicative character if one of these 
noises has many fewer time scale than others. In such a case, the limit cycle may 
be created if the most fast variable is coupled with more than two slow ones. However, 
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the case considered in Section H] shows that such dependence is not necessary to arrive 
at hmit cycle. The formal reason is that within adiabatic condition Tp <^ te both noise 
amplitude Qs{E) and generalized force J-'^^\E), determined with Eqs. ( l24l) and ( l25l) . 
enclose the squared strength E"^, but do not include the square S"^. 

The limit cycle is created if the fastest variations displays a principle variable, 
which is coupled with two different degrees of freedom or more. Indeed, at the relations 
Te <^ Tp, Ts of relaxation times considered in Section [3l the strength E evolves according 
to the stochastic law of motion ( |T7I) . Accounting this relation in the nonlinear terms of 
two last equations (fT^ arrives at dependencies of the noise amplitudes of the polarization 
P and the difference of level populations S on both variables S and P themselves. Due to 
gaussian nature of the noises their variances are additive values [39|B0]. so that effective 
noise amplitudes Qp and Qs of the principle variables are defined by Eqs. ( fT9l) . which 
include both squares 5*^ and P^. As a result, solutions of Eq.( l2Ti) become double- valued 
to be related to the limit cycle. 

This cycle appears physically as stochastic coherence, that has been observed 
both numerically [TH [31] and analytically [12] . Analogously to the doubly stochastic 
resonance [36], such a resonance may be organized if stochastic nonlinear system has two 
noises, but both of them must be multiplicative in nature. Moreover, the system under 
study acquires an intrinsic time scale related to multistable state only far off equilibrium 
statistical state. In opposite to the deterministic limit cycle which bounds a domain of 
unstable values of related variables, in our case stochastic variables evolve out off the 
limit cycle only, whereas in its interior the domain of forbidden values appears. 
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